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The angular, spectral and temporal features of the highest energy protons and accompanying them 
secondary neutrinos and synchrotron gamma-rays propagating through the intergalactic magnetic 
and radiation fields are studied using the analytical solutions of the Boltzmann transport equation 
obtained in the limit of the small-angle and continuous-energy-loss approximation. 
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I. INTRODUCTION 



Because of deflections in the interstellar and intergalac- 
tic magnetic fields, the information about the original 
directions of cosmic rays pointing to their production 
sites is lost. On the other hand, the isotropic flux of 
cosmic rays is contributed, most likely, by a large num- 
ber of galactic and extragalactic sources. These objects 
represent different source populations characterized by 
essentially different physical parameters - age, distance, 
energy budget, etc., as well as by different particle ac- 
celeration scenarios. This makes extremely difficult the 
identification of sources of cosmic rays based on the chem- 
ical composition and energy spectra of particles - two 
measurables characterizing the "soup" (isotropic flux of 
cosmic rays) cooked over cosmological timescales. For- 
tunately, at extremely high energies, E ~ 10 20 eV, the 
impact of galactic and extragalactic magnetic fields on 
the propagation of cosmic rays becomes less dramatic, 
which might result in large and small scale anisotropics of 
cosmic ray fluxes. Thus, depending on the strength and 
structure of the (highly unknown) intergalactic magnetic 
field (IGMF), the highest energy domain of cosmic rays 
may offer us a new astronomical discipline - " cosmic ray 
astronomy" . The extension of studies to energies 10 20 eV 
and beyond enhances the chances of localization of par- 
ticle accelerators for two reasons. With an increase of 
particle energy, the probability that a proton would pen- 
etrate through the intergalactic medium (IGM) without 
significant deflections in chaotic magnetic fields increases. 
Note that for IGMF much weaker than 1CT 9 G, the de- 
flection angle can be quite small also for lower energy 



'Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, 
Dublin 2, Ireland; Max-Planck-Institut fur Kernphysik, 
Saupferc hsckffiEg 1- D-6917 Heidelberg, Germany; Electronic 
address: 



Felix.Aharonian@mpi-hd.mpg.de 



t Moscow Institute of Engineering Physics, Kashirskoe sh. 31, 
Moscow, 115409 Russia; Max-Planck-Institut fur Kernphysik, 
Saupfe rcheckweg 1, D-6917 Heidelberg Germany ; Electronic ad- 
dress: 5tanislav~Kclncr@mDi-hd.mDg.dc] 

t Electronic address: Anton.Prosckin@mpi-hd.mpg.dc 



protons (9 cx B/E). However at energies significantly 
below 10 20 eV, the deflection in galactic magnetic fields 
becomes the dominant factor leading to the lost of in- 
formation about the original directions of particles (see, 
e.g., Ref. §). 

In the context of prospects of realization of "cosmic- 
ray astronomy" , there is a second independent factor 
which gives strong preference to energies 10 20 eV. Par- 
ticles of such high energies can arrive only from rela- 
tively nearby accelerators located within 100 Mpc (see, 
e.g., 0). This dramatically (by orders of magnitude) 
decreases the number of relevant sources of > 10 20 eV 
protons contributing to the observed cosmic ray flux, 
and correspondingly reduces the level of the diffuse back- 
ground, i.e. the (quasi) isotropic flux as a superposition 
of contributions by unresolved discrete sources. Formally, 
one cannot a priori exclude the possibility that the 10 20 
eV cosmic rays are contributed by a large number of weak 
sources which cannot be detected individually. Alterna- 
tively, the entire cosmic ray flux at such high energies can 
be dominated by contributions from a few sources, espe- 
cially given the tough requirements to the 10 20 eV pro- 
ton accelerators ||. This excludes, in particular, objects 
like ordinary galaxies, unless the galaxies provide high- 
est energy cosmic rays through transient events related 
to compact objects like Gamma Ray Bursts (GRB) Q. 

The propagation of cosmic rays in IGMF has been dis- 
cussed in a number of recent works (see, e.g., Refs. §- 
||). In these studies different magnetized environments 
have been assumed and different methods dealing with 
particle transport have be applied. Consequently, their 
conclusions are quite different, the principal reason being 
the different assumptions and approaches in the modeling 
of the IGMF. The main purpose of our work is the study 
of features related to the transport of particles, therefore 
we assume, following Ref. 0], purely turbulent and ho- 
mogeneous IGMF. This not only makes the calculations 
simple and more transparent (as long as it concerns the 
pure transport effects), but also seems to be a feasible 
realization for the large scale structure of IGMF. 

Whether we may identify the accelerators of extra- 
galactic cosmic rays using the highest energy protons is 
a question which largely depends on the strength of the 
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large scale IGMF. Even for the most favorable condi- 
tions for realization of the " proton astronomy" , the lat- 
ter will be relevant to the nearby Universe, the accessible 
sources being limited within a sphere of radius 100 Mpc. 
A different approach for localization of acceleration cites 
of 10 20 eV protons can be provided by observations of 
gamma-rays and neutrinos produced at interactions of 
these energetic particles with the 2.7 K CMBR pho- 
tons and magnetic fields in the proximity of the source, 
namely within a region of a size of order of 10 Mpc - 
sufficiently large for effective interaction of protons with 
2.7 K photons through photomeson process and, at the 
same time, still small for a significant deflection of pro- 
tons from their original directions. All short-lived par- 
ticles of these interactions, as well as the products of 
their decays (gamma-rays, neutrinos and electrons) are 
produced at small angles relative to the initial directions 
of parent protons. In an magnetized environment with 
B > 10 ~ 9 G the electrons with typical energy exceed- 
ing 10 19 eV are predominantly cooled via synchrotron 
radiation with production of high energy gamma-rays. 
The electrons emit synchrotron photons very quickly, be- 
fore any significant change of their direction in the sur- 
rounding chaotic magnetic field. Thus the synchrotron 
photons will move essentially in the initial direction of 
the parent protons. Since the protons, after they escape 
their production site (accelerator) , move radially, the ob- 
server will see an apparent compact (quasi-point like) 
gamma-ray source |}| |l^|, even though gamma-rays are 
produced in an extended region with angular size of or- 
der of ~ d/r ~ 5 (r/lOOMpc) -1 deg. Note that the same 
is true if protons escape the source anisotropically, but 
are moving within a narrow angular cone towards the 
observer. Otherwise, the observer will miss the source. 

The favorable range of IGMF for realization of this 
scenario is 10 -9 — 10 -8 G. In a stronger magnetic field, 
deflections of protons are significant even at the first sev- 
eral Mpc scales. Thus, because of the small interaction 
depth of undeviated protons, the point like source be- 
comes very weak. 

On the other hand, for IGMF much weaker than 10 _9 G 
electrons are cooled predominantly via inverse Compton 
scattering, thus the efficiency of synchrotron radiation 
drops dramatically. This scenario which involves a pair 
cascade in the 2.7 CMBR and Extragalactic Background 
Radiation (EBL), also leads to high energy gamma-rays. 
However, unless the field is much weaker than 10~ 12 G, 
the cascade electrons of relatively low (TeV) energies are 
thermalized, thus the cascade leads to the formation of 
giant halos |nj and in this way contribute to the dif- 
fuse extragalactic gamma-ray background radiation (see, 
e.g., Ref. p2| ) rather than to the formation of a discrete 
gamma-ray source (for a discussion of different regimes 
of formation of cascades initiated by interactions of high- 
est energy protons with 2.7 K CMBR, and their de- 
tectability from the direction of the cosmic ray source see 
Ref. jl3|). The detection of the cascade component as a 
point like or a slightly extended source of gamma-rays ini- 



tiated by interactions of ultrahigh energy protons (after 
they escape the accelerator) with 2.7 K CMBR is possi- 
ble in the case of extremely small IGMF, B < 10~ 15 G 
(see, e.g., Ref. @). 

The energy spectrum and flux of synchrotron radiation 
of secondary electrons from photomeson interactions of 
protons with 2.7 K CMBR have been studied in Ref. [jici|. 
The calculations have been limited by the first 10 Mpc 
range of propagation of protons, assuming that at this 
stage protons propagate radially without significant devi- 
ations, and the secondary electrons move along the same 
direction before they emit synchrotron photons. While 
this approximation gives a correct estimate of the flux, it 
does not specify the angle within which the radiation is 
confined. This approach ignores also the non- negligible 
tails of distribution of synchrotron radiation formed at 
the later stages of propagation and interactions of pro- 
tons. 

In the case of quasi-continuous operation of an extra- 
galactic accelerator of protons over timescales exceeding 
the typical delay time due to the deflection in the mag- 
netic field, the energy and angular distributions of pro- 
tons, as well as accompanying photons and electrons, can 
be accurately described by the steady-state solutions of 
the transport equations. Generally, this is the case of a 
continuous proton accelerator of age T > 10 6 yr. In the 
case of shorter activity of the source (an "impulsive ac- 
celerator") or solitary events like gamma-ray bursts, rel- 
atively simple analytical solutions of the arrival time dis- 
tributions of protons, gamma-rays and neutrinos can be 
obtained within an approximation when the energy losses 
of protons are ignored. We consider the cases of "contin- 
uous" and "impulsive" proton accelerations in Sections 
and IV, respectively. 



II. 



STEADY STATE DISTRIBUTION 
FUNCTIONS 



The realization of the small-angle multiple scattering 
considerably simplifies the description of propagation of 
protons through a scattering medium. In particular, in 
the small-angle approximation the term v |£ of the Boltz- 
mann transport equation can be presented in a form al- 
lowing analytical derivation of the steady state solution. 
Because of smallness of the single scattering angle one can 
write the elastic collision integral in the Fokker-Planck 
approximation. To expand the distribution function into 
series in terms of the single scattering angle one should 
have a smooth function of this angle. This condition is 
satisfied if one neglects unscattered part of the distri- 
bution function that has very sharp angle dependency. 
Such an approximation is justified in the case of multiple 
scattering. 

The approach provides solutions that can be applied 
to the various cases which, independent of the details 
of the scattering medium, are characterized only by the 
average scattering angle per unit length (9 2 s ). The scat- 
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tering process depends on the particle energy, i.e. (9 2 ) 
is a function of energy. During the propagation through 
the medium between two scattering centers, the energy 
of particles is gradually decreased due to different dissi- 
pative processes. If the change of energy in each action 
of interaction is considerably smaller than the initial en- 
ergy, one can use the continuous energy loss approxima- 
tion. It should be noted that in the approach described 
here the processes responsible for the scattering and the 
energy loss of particle are not required to be the same. 
The particle scattering could have elastic character and 
do not cause energy losses. On the other hand, the ef- 
fect of deflection of particles from their original direction 
due to the processes responsible for energy losses might 
be negligibly small. This is the case of the problem con- 
sidered below. One can safely ignore the change of the 
direction of primary particles as well as the production 
angles (6 ~ 1/7) of the secondary products (gamma-rays, 
electrons, neutrinos) due to all relevant processes includ- 
ing photo-meson and pair production, inverse Compton 
scattering, synchrotron radiation. 

The aim of this section is to derive distribution func- 
tions for protons and accompanying them secondary par- 
ticles propagating through the galactic and extragalactic 
magnetic fields for a spherically symmetric point source 
of protons. However, it is technically more convenient to 
consider first a source emitting protons in a given (fixed) 
direction. In this case we have a preferential direction 
along the infinitely narrow beam emitted by the source. 
Let us choose z-axis along this direction. Because of the 
scattering, particles deviate from the initial course. To 
define the deviation we introduce angles 9 X and 9 y be- 
tween the direction of propagation n and the coordinate 
planes YOZ and XOZ, respectively. If n is close to z- 
axis, the angles 9 X and 9 y are small and can be treated 
as components of two-dimensional vector 9 in the XOY 
plane, where the absolute value of 9 corresponds to the 
deflection angle between n and the z-axis. Then we can 
write n a (9 x ,9 y , 1 - f ) = (0, 1 - £). 

The retention of the second-order term 9 2 /2 in the ex- 
pansion of n z allows us to take into account the effects 
relating to the elongation of the path like delay time, but 
does not give any considerable contribution to steady- 
state solution. Therefore we divide the problem into two 
sub-problems. In the first part of the paper we solve the 
steady-state equation that takes into account the energy 
losses but ignores the elongation of particle trajectories. 
The results of these calculations are relevant to the " con- 
tinuous" source of protons and describe the energy and 
angular distributions of protons and accompanying neu- 
trinos and synchrotron radiation of secondary electrons 
produced during the propagation of protons. In the sec- 
ond part of the paper we calculate the distributions of 
arrival times of protons, neutrinos and gamma-rays in 
the case of an "impulsive" source. In this case the ar- 
rival time delays directly depend on the elongation of 
trajectory. The time-dependent solutions for distribu- 
tion functions presented in Section IV are limited by the 



approximation in which the energy losses of protons are 
neglected. 

Thus, to derive the steady-state solution of the trans- 
port equation we assume n « (9 x ,9 y ,l) = (6,1) . Let 
us denote by p = (x, y) the perpendicular displacement 
in the plane XOY. For a point source characterized by 
a monoenergetic and infinitely narrow beam of protons 
emitted along the z-axis we obtain the equation for Green 
function G(r, 6, E) of the Boltzmann steady-state trans- 
port equation in the approximations of a small-angle mul- 
tiple scattering and continuous energy losses: 



d_ 

dz 







dp 







e)G{r,9,E,E a ) 



4 89 2 dE 

= h(z)S(p)6(e)6(E-E ). (1) 



Here we take into account that the particles are ultrarel- 
ativistic \v\ = c. The solution of Eq. (nl) is obtained in 



Rcf. [15 for the propagation of charged particles passing 
through a layer of matter. The features of this solution 
are comprehensively discussed in Ref . Jl(| . Using the no- 
tations introduced in Ref. [M, the Green function can 
be written in the form: 



G(r,9,E,E Q ) 
A lP 2 - 



x exp 



_ S(S(E,E )-z) 

ce(E)ir 2 A 
2A 2 9p + A 3 9 2 



(2) 



where S is the traveled distance that is uniquely related 
to the energy loss rate e(E) — \dE/dz\: 



S(E,E ) 



dE' 



and 



A = A X A 3 - A\ 



(3) 



(4) 



The ^-function in Eq. (Q) points to the fact that we ne- 
glect the elongation of trajectory so the traveled distance 
is equal to z as if particles propagate strictly along z-axis. 
Taking the relation between energy and z into account, 
Ai can be written in the following form: 

z 

ME ,z) = J (9 2 )(z')(z-zy- 1 dz'. (5) 


It is easy to recognize the physical meanings of the co- 
efficients A\, A2 and ^3; A\ is the mean square deflection 
angle, A3 is the mean square displacement, and A2 is the 
mean value of 9p at the distance z: 



Ai = (9 2 ) z , A 2 = (9p) z , A 3 = (p 2 



(6) 



For the treatment of the case of spherically symmetric 
point source of protons, let us rewrite Green function in 



4 



the form which is independent of choice of the coordinate 
system. After the replacements 



where no is the direction of the emission, n is the direc- 
tion of particle motion at the point r, we find 



e 



n - n , 



rn , 



r, 



(7) 



5(S(E,E )-r) ( A 1 (r-rn ) 2 -2A 2 (r-rn )(n-n a )+A 3 (n-n ) 2 \ 



(8) 



Performing integration over all directions of the vector 
no by the saddle point method (see Appendix |a|), we 
find ' 



where 



D = Ax 



.Ao 



A 



(10) 



Since we have spherically symmetric distribution, the 
Green function depends only on the distance r from the 
source and 6 which is the angle between the radius- vector 
from the source to the observation point and the move- 
ment direction at this point. 

We assume that the spherically symmetric source in- 
jects protons into the intergalactic medium with a con- 
stant rate: 



Q p (r,E) = cJ p (E)S(r). 
The substitution of this expression into 



(11) 



f{r,9,E) = J Q(r ,E ) 
xG sph (r - r ,e,E,Eo)dr dEo 



(12) 



gives 



f P (r,0,E) 



Jp(E ) 
e(E) J 7rr 2 D 

E 



exp 



D 



xS(S(E,E )-r)dE , (13) 



where D can be written as 



D(E ,r) 



(0 2 s)(r r y 2 dr>. 



(14) 



For the given energy and spatial distribution of pro- 
tons we can calculate number of secondary particles from 
the decays of 7r-mesons that are produced at interactions 
between protons and 2.7 K CMBR photons. To obtain 
the energy distributions of the secondary products - pho- 
tons, electrons and neutrinos, we use the approximation 



proposed in Ref. [[l7| . The energy of protons is ultrarela- 
tivistic so we can assume that secondary particles initially 
move in the same direction as protons. The distribution 
of second particles can be presented in the form 



Q(r,6,E) = Q(f p (r,0,E p )) 



(15) 



where Q denotes an integral operator. For example, for 
the energy distribution of protons J P (E), the energy dis- 
tribution of photons produced in photomeson interac- 
tions is 



f(E-y) = Qj(j p (E p )), 



(16) 



where 



Q 1 {J p {E p )) = J J p (E p )f ph (e)w(E 7 ,E p ,e)dE p de . (17) 

Here f p h is the distribution function of CMBR photons, 
w is the differential interaction rate of the 737 interactions, 
namely, the Bethe- Heftier pair production or photomeson 
production (see Ref. Q). Since we are interested in 
the distribution of ultrarelativistic electrons that weakly 
deviate in the magnetic field, we can apply the Green 
function given by Eq. (|) to the source function given 
by Eq. (|l5|). Note that Q acts only on variable E p , 
therefore we can change the order of integration. Tedious 
calculations (see Appendix ph yield: 



fe{r,9,E e ) 



1 



ce e (E e ) 



dE e0 Q c 



6 p( e p) 



dE, 



J p (E p0 )^{-D^D p 

X ^ n(D e + D p ) 6{S ?) 



(18) 



Here S is sum of the distances traveled by proton to the 
point of interaction with CMRB and traveled by electron 
from the point of production to the point r: 



S — S p (E p , E p o) + S e {E e , E e o). 



(19) 



The angular distribution of electrons in Eq. (|l^) is char- 
acterized by 



D P 



A e l 



(20) 
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where A e i = A e i(E e , E e o) have the same meaning as in 
Eq. (§), and 



D 



(21) 



where r = S p (E p ,E p0 ). 

The main channel of production of gamma rays by HE 
electrons is s ynchrotron radiation. Applying the modi- 
fied Eq. (Dl) for the spectrum of synchrotron radiation 
in chaotic magnetic field to the distribution of electrons 
given by Eq. (|l8|), we find the angular and spatial distri- 
butions of gamma ray. Let us describe the procedure as 
in Eq. © by 



,(r,9,E 7 ) = Q s (f e (r,9,E e )). 



(22) 



In the general case the distribution function of gamma 
rays that are characterized by the source function 
Q(r, n, E) and propagate through the absorbing medium 
with extinction coefficient k(E), is 



fr 1 (r,n,E) = ~f Q(r - nT,n,E)e- kT dT, 



(23) 



we obtain: 



f 1 {r,6,E 1 ) = Q s 



c 2 e e (E e ) 



dE e oQ e 



x J dE p of p Ti(l - ^\ rjf^ - sin( 



, (26) 



where 



f P = 



Jp{Epo) 



exp 



D c + D 



, (27) 



r V V r , 



r\ is Heaviside function. The angle between n and r — nr 
is 



'(S/r) 2 - sin 2 6 



i\ = arccos 



S/r 



(28) 



III. THE SPECTRAL AND ANGULAR 
DISTRIBUTIONS OF PROTONS, PHOTONS 
AND NEUTRINOS 

A. Protons 



where n is the direction of movement at the point r that 
coincides with the direction of emitting as the propaga- 
tion of gamma rays is rectilinear. Since emitting electrons 
have ultrarclativistic energy we assume that the direction 
of radiation coincide with the electron direction, there- 
fore the distribution function of synchrotron gamma rays 
for the distribution of electrons given by Eq. (18) is 



/ 7 (r, 9, = ~ J Q s (M\r - nr\A, E e )) e~ k ^ T dT, 
o 

(24) 

where Oi is the angle between n and r — nr, and 9 is 
the angle between n and r. It is convenient to perform 
the integration over r using delta- function in Eq. (18). 
For that we should change the order of integration so the 
integration over r becomes internal. Using the features 
of the delta-function we find 



*(|r 



TIT 



S) 



\r — nr 



Sr 



(f) " 



sm 



2 

1 i=i 



5(i 



(25) 

where Tj.,2 = ^(cos^ =p \J ) 2 — sm2 Since the term 
corresponding to r 2 does not contribute to large angles 
8i in exponent (see Eq. ([l8])) we keep only the term cor- 
responding to T\ . After performing relevant calculations 



Transport of protons substantially depends on the spa- 
tial distribution of magnetic fields. The assumption of 
chaotically oriented magnetic cells is usually used for es- 
timates of the influence of IGMF on cosmic-ray propa- 
gation (see, e.g., Ref. fl8| ). The spectral analysis of the 
correlation function of the magnetic field fluctuations |lJJ] 
provides a more appropriate and accurate treatment of 
the problem. We use this approach for derivation of the 
mean square deflection angle. 

UHE protons propagate large distances in IGMF with- 
out considerable deflections. Indeed, the evaluation of 
the deflection angle 89 ~ X/r g on the correlation length 
A is 



69 ~ 9x10" 



A 



IMpc 



B 



10" 9 G 



10 20 eV 
E~ 



rad, (29) 



where r g = E/eB is the gyroradius of the ultrarelativistic 
particle. Therefore change of direction of ultrahigh en- 
ergy protons is small on the scale A ~ 1 Mpc. The proton 
energy can be assumed constant for this scale. Then the 
proton motion in the magnetic field is described by the 
equation 



ec 
E 



[v x B(r)} . 



(30) 



For ultrarelativistic particles v — cn, where n is a unit 
vector. Rewriting the change in velocity over the time 
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At in the form Av — c9, we find 



t+At 



6= J ^[nxB(r(t))]dt. 



(31) 



Since the deflection angle is small, the trajectory of par- 
ticle can be considered rectilinear in integration. 

Now one should make an assumption about the statis- 
tical properties of the magnetic fields. Here we assume 
that IGMF is a statistically isotropic and homogeneous. 
While (6) = (since in this case (B) = 0), the mean 
square deflection is 

( d2 ) = (^Y f([nxB 1 ][nxB 2 })dt 1 dt 2 



= (jjj) (<W - n a n a ) j (B la B 2 p) dz x dz 2 , (32) 

where B\ f2 = B(r(ti j2 )). Here we switch to integration 
over coordinates of particle, directing z-axis along n. 

Eq. (|3^) includes the correlation function of the mag- 
netic field 



K aP iTx-T 2 ) = (B a ( ri )Bp(r 2 )} 



(33) 



It depends only on the difference (ri — r 2 ) because we 
assume statistical homogeneity of the magnetic field. The 
mean square of magnetic field is determined as (B 2 ) = 
K aa (0) = const. 

To turn to the spectral description, K a p should be 
written as a Fourier integral: 



K a0 (ri -r 2 )= I K aj3 {k)e 



a ik{r 1 -r 2 ) 



d 3 k 

J2nf 



(34) 



Since divB — 0, then K a p should satisfy the conditions 
dK aj3 /dxi a = , dK a/ 3/dx 2l3 = . (35) 

For function K a p these conditions take on form K a pk a — 
0, K a fjkp = 0. Therefore, if there is no preferential di- 
rection in space, K a p has the following structure: 



$(fc 2 ) (B 2 ) . (36) 



Here constant factor (B 2 ) is introduced such that $(fc 2 ) 
meets the normalization condition: 



$(fc 2 



d A k 
(2nf 



OO 

{ <S>(k 2 )k 2 dk = 1. (37) 

l-K 1 J 



It is convenient to change variables z\ = z + Q/2, 
z 2 = z — (/2 in the integral in Eq. (|32|). Assuming that 
the traveled distance Az is much greater than a char- 
acteristic scale on which the correlation function tends 
to zero, one can extend the limits of integration over £ 



to infinity. Meanwhile, the traveled distance should be 
smaller than the distance on which the proton loses its 
energy appreciably. The integrand depends only on (, 
therefore the integration over dz gives the length of inte- 
gration interval Az. The mean square deflection angle is 
proportional to the traveled distance, so the mean square 
deflection per unit length is 



Jf a/J (0,0,C)dC. (38) 



(J;) (S a p - n a n p ) 



Integration of K a p written in the form of Eq. (|34|) over 
d( gives 2TrS(k z ) that allows us to calculate the integral 
over dk z . Then, the integral remains over the components 
of k perpendicular to z-axis. Using Eq. (dw, we obtain 



1 / e \ 2 



d 2 kj_ , , 
2^ - / (2^ $( ^ 



(B 



The derived result can be written in the form 

,„o, 7T / e n 2 



2 \E 



(B 2 )A. 



(39) 



(40) 



where 
A = 



Taking into account Eq. ( p7|) , the factor A can be treated 
as the mean value of A = (fc _1 ). 

The calculation of (9 2 ) requires the spectral energy dis- 
tribution of magnetic field. To obtain final form of (6 2 s ) 
we assume a power-law spectrum 



$(fc 2 )fc 2 



> k<k o 



(42) 



where fco is an absolute value of the wave vector corre- 
sponding to the maximal scale of correlation A: fco = 
2tt/X. It gives 



4a/3 



(l) ; 



(B 2 )X. 



(43) 



Taking into account the turbulent character of IGMF 
that has (B) = 0, we take a = 5/3 which corresponds to 
the Kolmogorov turbulence. The choice of the parameter 
/3 is, to a certain extent, arbitrary. Here we assume j3 — 1 
which leads to a simple expression for the mean square 
deflection angle per unit length: 



(44) 



Because of uncertainties related to the spectrum of 
IGMF, the numerical factor in Eq. ( [44| ) is somewhat dif- 
ferent from the coefficients used in other papers (see, e.g., 
Ref. H). 
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Significant uncertainty in calculations of (9 2 ) is related 
to the absolute value of the correlation length A. It is ex- 
pected to be between 100 kpc and 1 Mpc, i.e. comparable 
to the the characteristic distances between galaxies. In 
the subsequent calculations we normalize the correlation 
length to A = 1 Mpc, but the presented results can be 
easily recalculated for any A. 

Since the propagation of protons in IGMF can be 
treated as a set of large number of small chaotic deflec- 
tions, the problem can be reduced to the diffusion in an- 
gle. The diffusion coefficient D{r, E ) given by Eq. (|i~4|)) 
contains information about the energy loss and influence 
of IGMF on propagation, and gives angular distribution 
of protons at the given point. Since there is a unique cor- 
respondence between the energy and r (see Eq. (^)) we 
can rewrite Eq.([wj) in terms of energy and energy losses 
per unit length. Substituting Eq.([44|) into Eq. (|l4|), we 
obtain 



D(E,E ) 




dE' 



(45) 



where 



E 



dE' 



- 2 A 



(46) 




FIG. 1: The mean free path of protons in the intergalactic 
medium due to interactions with photons of CMBR. 



The function E/e(E) based on results of Ref. |17| and 
implying the mean free path of protons in the intergalac- 
tic medium due to the Bethe-Heitler pair-production and 
photomeson processes at interactions with CMBR, is 
shown in Fig. [|. 

Note that for many scenarios described by Eq. (|l|) the 
same process is responsible for both the angular scatter- 
ing and the energy losses. But in the case of propagation 
of protons in the intergalactic medium we deal with two 
different processes: while the interactions with CMBR 
lead to energy losses, the angular deflections are caused 
by multiple scattering on magnetic inhomogeneities. 

The influence of energy losses on the angular distribu- 
tion of protons can be traced in Fig. [2| where mean 
deflection angle of protons with observed energies Ef 
is shown as function of traveled distance r. As it is 
seen from Fig. [j], protons with energy smaller than 
E = 6 x 10 19 eV do not suffer noticeable energy loses 
over the distances ~ 100 Mpc. In this case the diffu- 
sion in angular space can be treated as a homogeneous 
random walk that brings us to the dependence of the 
mean deflection angle on the travel distance oc r 1 / 2 . For 
protons with initial energy higher than the threshold of 
photomeson production, the energy of protons gradually 
decreases which leads to deviation from this simple de- 
pendence. In particular, for the given observed (final) 
energy Ef, this effect implies higher original energies, 
and consequently smaller deflection angles at the initial 
parts of propagation. This results in a weaker increase of 
the mean deflection angle with the traveled distance in 
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FIG. 2: The mean deflection angle of protons for the fixed 
observed energy Ef over the distance r. The numbers at the 
curves indicate the energies which proton had at the distance 
r from the observer. 



comparison with loss-free case. This effect is clearly seen 
from analytical expressions, which is possible to obtain in 
the case of constant energy loss rate 1 4^ I -= = b = const: 



E 2 



C 2 - 2C + 2(1 



(47) 



\C=2br 
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Expanding this expression into series in terms of powers 
of r we obtain: 



E) V3 



br^ 
6 



(48) 



The first term does not depend on the value of b and 
thereby describes the loss-free propagation. The next 
term takes into account the energy losses and makes the 
dependence on the distance r weaker. While Eq. ( |47| ) 
approximately describes the behavior of mean deflection 
angle for the final energy Ef > 10 21 eV, the first term 
of Eq. (48|) describes the case of Ef < 6 x 10 19 eV. In 
Fig. U the mean deflection angle of protons is given for 
IGMF B = I nG. Since the dependence of the average 
deflection angle on the magnetic field is linear, it is easy 
to produce plots for other magnetic fields. 

In order to indicate the evolution of the energy of pro- 
tons during their propagation through the 2.7 K CMBR, 
in Fig. ^ we indicate at the corresponding curves, cal- 
culated for the fixed observed (final) energies of protons, 
the energies which protons had at different distances from 
the observer. For the fixed observed energies exceeding 
the threshold of photomeson production, the calculated 
initial energies grow dramatically with the increase of the 
distance, especially for > 100 Mpc. Therefore any deficit 
of protons of such high energies in the initial spectrum 
would results in the cutoff in observed spectrum at the 
corresponding energies. 

The energy distributions of protons at distances 100 
Mpc and 300 Mpc are shown in Fig. || for the initial dif- 
ferential energy spectrum J P (E) = JqE~ 2 cxp(— E/E ). 
The total luminosity of the source in CRs with energy 
above 10 9 eV is taken L = I0 44 erg/s. The upper dashed 
lines correspond to the case when protons propagate in 
empty space; flux is determined by the geometrical fac- 
tor 1/r 2 . Solid line presents the case when the deflections 
in the magnetic field are ignored. Comparison of these 
two curves reveals two features: a bump and a sagging 
at lower energies. Both features become more prominent 
with increasing of the distance. The bump preceding the 
cutoff appears due to strong growth of energy losses at 
the threshold of photomeson production (see Fig. [j]) that 
makes particles to be accumulated in this energy region; 
the sagging is a consequence of the energy losses due to 
the electron-positron pair production (see, e.g., Ref. [pp[). 

The approximation of continuous energy losses takes 
into consideration the mean energy losses. In general it 
provides an acceptable accuracy but some features con- 
nected with stochastic properties of interactions should 
be taken into account for precise description of the spec- 
trum in the cutoff region. The fluctuations in the energy 
losses do have an impact on the form of the bump and the 
cutoff in the observed spectrum of protons. It results, in 
particular, in a smoother cutoff and a broader and lower- 
amplitude bump |23| compared to the results calculated 
within the continuous energy losses approximation. 

The impact of the magnetic field leads to strong de- 
pendence of the energy distribution on the solid angle 



within which the particles are detected. As it is seen in 
Fig. |[ the flux of protons at highest energies is concen- 
trated along the direction to the source; the protons of 
lower energies are scattered over large angles. 



B. Electrons 

The secondary gamma-rays and neutrinos are tracers 
of propagation of protons in the intergalactic medium. 
The first generation gamma-rays from photomeson pro- 
cesses are produced at extremely high energies E > 
I0 19 eV. They are effectively absorbed due to interac- 
tions with the photons of CMBR and the Cosmic Radio 
Background (CRB) over distance ~ 1 Mpc. Because of 
the threshold effects, at energies below 10 14 eV the ef- 
ficiency of interactions with CMBR dramatically drops, 
but gamma-rays continue to interact with the infrared 
and optical photons of the Extragalactic Background 
Light (EBL). At these energies the mean three path 
of gamma-rays increases sharply achieving, ~100 Mpc 
at £ 7 -10 TeV, and -I Gpc at E < 200 GeV (see, 
e.g., Ref. |pi[). Thus, as long as we are interested in 
gamma-rays from the sources of highest energy cosmic 
rays, the energy of gamma-rays should not significantly 
exceed I TeV. In this energy band gamma-rays are pro- 
duced through the electromagnetic cascade initiated by 
the products of decays of short-lived mesons from the 
photomeson interactions and, partly, by electrons from 
the Bethe-Heitler pair-production process. For the devel- 
opment of an effective cascade the magnetic field should 
be smaller than I0~ 10 G. Even so, the observer can see 
the cascade gamma-rays in the direction of the source 
only in the case of extremely small IGMF, B < I0 _15 G. 
A collimated beam of gamma-rays of GeV-TeV ener- 
gies is expected in the case of magnetized intergalactic 
medium with B > I0 _9 G. These gamma-rays are pro- 
duced through the synchrotron radiation of E > 10 19 eV 
electrons. 

Due to very large Lorentz factor of particles, we can 
assume that secondary products from all interactions un- 
der consideration propagate strictly in the direction of 
the parent particle. Therefore observed angular distribu- 
tion of gamma rays depends on the influence of IGMF 
on electrons that produce these gamma rays. To observe 
the UHECR source in gamma rays it is necessary that 
producing electrons are only slightly deflected in IGMF. 

The almost rectilinear part of the path of electrons is 
much smaller than distances traveled by protons and is 
comparable to the typical correlation length, A ~ I Mpc. 
So the scattering of electrons takes place in almost homo- 
geneous magnetic field. But since direction of magnetic 
field have a random character the scattering occurs in 
random directions. In case of electrons one can apply 
the formalism of the multiple scattering to the random 
single scattering. Indeed, as have been noted the distri- 
bution function should be smooth function of angle to 
write the elastic collision integral in the Fokker-Planck 
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FIG. 3: Energy flux distribution of protons observed within different angles for the source at the distance r = 100 Mpc (left 
panel) and r = 300 Mpc (right panel). The initial spectrum of protons is assumed power-law with an exponential cutoff at 
E = 3 x 10 20 eV, the IGMF is 1 nG. 



approximation. If all particles are scattered, as in the 
case of electrons, the distribution function does not in- 
clude a part with sharp angular distribution correspond- 
ing to non-scattered particles. 

To obtain the mean square deflection angle per unit 
length we use expression for deflection angle of ultrarel- 
ativistic electron traveled the path on which its energy 
has changed from the initial energy E' to the hnal energy 
E: 



E' 



-dE' 1 



e e r g 



(49) 



where e e is rate of energy losses due to synchrotron radi- 
ation and r g is gyroradius. Taking into account random 
field orientations we find 

3 (^ c2 ) 4 









2 1 \ 




1) 




h 3^) 



(50) 



After substitution this equation into Eq. we find the 
coefficients A e i of the diffusion coefficient D e in simple 
analytical forms: 



A 1 



A 2 = 



5cv/3 1 
180^4 



(7^ 2 + 14^ + 9)(i - e) 2 



(51) 



q/3 2 1 
180 E 5 



(19£ 2 + 22£ + 9)(l-£) 3 , (52) 



As = j^-^ (12C 2 + 10£ + 3)(1 - £) 4 , (53) 



where 



6 (mc 2 ) 
a = — 



Eq. 



9 (mc 2 ) 4 

5 e 2 ' " ~ 4 e 4 B 2 ' 
can be written in the form 



e « 0.008' 



B nG El 



(54) 



(55) 



where B n c is the magnetic field in units of nanoGauss 
(nG), E 2 q is final energy in units of 10 20 eV, £ = E/E'. 
Here the random orientations of the field are taken into 
account. This expression allows us to estimate the 
threshold of isotropization. Indeed, if the electron loses 
considerable part of its energy, then £ -C 1 and deflection 
angle mostly depends on the final energy. The deflection 
angle becomes quite large (~ 1 radian) in the magnetic 
of field 1 nG when final energy is E ps 2 x 10 18 eV. For 
greater magnetic field the threshold of isotropization is 
shifted to the range of lower energies. It should be noted 
that for magnetic fields 1 — 100 nG this threshold appears 
in the energy region where the energy losses due to syn- 
chrotron radiation dominate over the inverse Compton 
scattering (see Fig. ^). It means that inverse Compton 
scattering can be neglected for electrons under consider- 
ation. 

Let us estimate the energy of gamma rays produced 
by electrons with energy exceeding the threshold of 
isotropization. Using modified Eq. (Dl) for energy dis- 
tribution of synchrotron radiation in chaotic magnetic 
fields, we find the energies of electrons that produce syn- 
chrotron gamma rays with energy Ey. 



E e = 1.23 x 10 1 



xB. 



nG 



(56) 
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FIG. 4: Energy loss rates of electrons due to inverse Comp- 
ton scattering on CMBR photons (solid line) and synchrotron 
radiation in random magnetic field for B =1 nG, 10 nG, and 
100 nG. For electrons of energy E > 10 19 eV the inverse Comp- 
ton scattering on the radiowaves of CRB becomes comparable 
or even can exceed the contribution of the Compton scatter- 
ing on CMBR, however for IGMF B > 1 nG the synchrotron 
radiation remains the main cooling channel. 



Here E e and E 7 are given in units of eV, x is the dimen- 
sionless argument of distribution function Eq. (D5). The 
latter has a maximum at x 



0.2291 and exponentially 
decreases for large x (see Eq. (D7)). To make sure that 
observed gamma rays are produced by electrons with en- 
ergies greater threshold of isotropization we should con- 
sider gamma rays with energies E 1 > 10 9 eV. Indeed, 
electrons with energies corresponding to i > 10 in the 
Eq. ( p^ ) give exponentially small contribution into radi- 
ation of gamma rays of the given energy. Therefore, as- 
suming x = 10, we find that the contribution of electrons 
with energies below threshold of isotropization E e < 10 18 
eV into radiation of = 10 9 eV gamma rays is insignifi- 
cant. According to Eq. ( |55| ) the product BE 2 is constant 
for the isotropization threshold. Since the same combi- 
nation enters in Eq. ( |56| ) the minimal energy of gamma 
rays produced by the electrons under consideration does 
not depend on the magnetic field. 

At interactions of protons with the intergalactic radi- 
ation fields the ultrahigh energy electrons are produced 
via two channels: pair production and photomeson pro- 
duction processes. In the pair production process only a 
small (< 2m e /m p ) fraction of proton energy is converted 
to the secondary electrons. For the magnetic field of or- 
der of nG or larger, the energies of these electrons appear 
below the threshold of isotropization, thus they do not 
contribute to the gamma-ray emission emitted towards 
the observer. The photomeson processes lead to several 
non-stable secondary particles, such as tt, rj, K mesons, 



FIG. 5: Number of electrons of energy E e located inside a 
sphere of the radius r. 



which decay into high energy gamma rays, neutrinos and 
electrons. The electrons from the decays of these mesons 
are produced with energies flTj exceeding the isotropiza- 
tion threshold. 

In addition, a significant fraction of electrons is cre- 
ated at interactions of the first generation ("photome- 
son") gamma-rays with photons of CMBR and CRB. For 
the model of CRB suggested by Q , the mean free path 
of gamma-rays of E > 10 19 eV is determined by the in- 
teractions with MHz radiowaves; it is of order of several 
Mpc. Here we neglect by the interaction length assuming 
that gamma rays interact with CRB immediately after 
their creation. In this case the particle get additional de- 
flection since it is treated as electron all along. It results 
in broader angular distribution of observed gamma ray 
in comparison with exact consideration. The interaction 
of gamma rays with CRB photons of energy eji occurs 
in the regime Cfli? 7 /mgC 4 3> 1. It means that the most 
of the energy is converted to one of the two electrons. 
The energy of gamma rays is higher than the energy of 
electrons produced in the decays of mesons (see, e.g., 
Ref. [[l7]). Therefore electrons created by pair produc- 
tion process are more energetic than electrons generated 
in the decays of nonstable products of photomeson pro- 
cesses. Consequently, the pair- produced electrons result 
in higher flux of synchrotron radiation than the direct 
ones from the meson decays. 



C. Gamma rays and neutrinos 

The apparent angular size of the synchrotron gamma- 
ray source depends on the linear size of the emitter itself 
and the deflection angles of the parent electrons. Both 
are denned by spatial and angular distributions of elec- 
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trons, respectively. In the case of spherically symmet- 
ric source and small deflection angles of electrons 6d e f, 
the source located at the distance r with the gamma- 
ray emission region of radius d, has an angular size 
•dobs ~ 2^9def- The case of isotropically emitted gamma- 
ray source corresponds to Odef ~ 1- The linear size of 
the gamma-ray emitter can be evaluated from Fig. |^, 
where is shown the number of electrons of energy E e lo- 
cated inside the sphere of radius r. The saturation that 
takes place at large distances shows the absence of elec- 
trons in this region. One can see from Fig. [5] that the 
size of the sphere, where the electrons are located, de- 
creases while the energy increases. This is explained by 
the fact that protons producing electrons of such high en- 
ergies disappear due to energy losses. For energies below 
E e — 10 19 eV the electrons are not located in a definite re- 
gion. The electrons with energy below the thermalization 
threshold form an extended halo. These electrons have 
energies at which the inverse Compton scattering losses 
dominates over the energy losses due to the synchrotron 
radiation. They initiate electromagnetic cascades in the 
CMBR and EBL photon fields that eventually results in 
a very extended GeV-TeV gamma ray emission. 

The spectral energy distributions (SED) of gamma- 
rays, E 2 F(E), received within different angles are pre- 
sented in Fig. [| The fluxes are calculated for the same 
initial proton energy distribution used in Fig. [| Three 
series of curves for each of two distances (left and right 
panels) correspond to different cutoff energies in the ini- 
tial proton spectrum. One can see that the cutoff en- 
ergy has significant impact on the flux of gamma rays; it 
increases the flux, shifts the maximum of SED towards 
higher energies, and makes narrower the angular distri- 
butions. These features have a simple explanation. The 
increase of the cutoff energy provides more secondary 
electrons and extends the spectrum of electrons to more 
energetic region. The latter leads to smaller deflections. 
It is interesting to note that although the angular dis- 
tribution of gamma rays is composed of deflections of 
both protons and electrons, their angular distribution is 
more narrow compared to the angular distribution of pro- 
tons (see, Fig. ||). This is explained by the fact that the 
main portion of gamma rays is produced in regions close 
to the source by the highest energy protons which did 
not suffer significant energy losses (see, Fig. ||), while 
the multiple scattering in IGMF contributes to the for- 
mation of the angular distribution of protons over the 
entire path from the source to the observer. Since the 
angular size of the gamma-ray source is determined by 
the geometrical factor d/r, the distribution of gamma- 
rays from a source at the distance r — 300 Mpc is nar- 
rower than from an identical source located at the dis- 
tance r = 100 Mpc. It is remarkable that at very high 
energies the source becomes point-like. In particular, at 
energies above E 1 sw 5 x 10 11 eV, the observer will see 
the gamma-ray source located at the distance of 100 Mpc 
within an angle smaller than o b s = 0.1°. 

Fig. [7] shows the impact of the IGMF strength on the 



flux distribution of gamma rays. The increase of the mag- 
netic field leads to the shift of the maximum of SED to 
higher energies. In accordance with Eq. (E^), the shift of 
the synchrotron peak is proportional to the strength of 
the magnetic field since the energy distribution of elec- 
trons does not depend on the magnetic field. Finally 
note that the increase of the magnetic field implies strong 
deflections which leads to the reduction of the flux and 
widening of the angular distribution of gamma-rays. 

For the sources located beyond 100 Mpc, TeV gamma- 
rays interact effectively with optical and infrared pho- 
tons of the Extragalactic Background Light (EBL). The 
energy- dependent absorption of gamma-rays is character- 
ized by the optical depth r 77 which depends on the EBL 
flux and is proportional to the distance to the source. 
Unfortunately the EBL flux contains quite large uncer- 
tainties, especially at the mid and far IR wavelengths 
which are most relevant to the gamma-ray energy band 
and the source distances discussed in this paper. The im- 
pact of these uncertainties on the intergalactic absorption 
of gamma-rays is discussed in Ref . |2lJ . Even for the min- 
imum EBL flux at infrared wavelengths, the absorption 
of TeV gamma-rays from sources beyond 100 Mpc can be 
significant; at multi-TeV energies the optical depth t 77 
exceed 1 . Therefore the curves in Figs ^ , [?] and |s| should 
be corrected by multiplying the unabsorbed fluxes to the 
factor exp(— r 77 ). 

The decay of nonstable products of photomeson pro- 
cesses leads to the appearance of extremely high en- 
ergy electrons (positrons) and neutrinos (antineutrinos) . 
Since the magnetic field does not have an impact on neu- 
trinos, the angular distribution of neutrinos is determined 
only by the deflection of protons. This leads to more nar- 
row angular distributions of neutrinos compared not only 
to the distributions of protons (for the same reason de- 
scribed above for gama-rays) but also compared to the 
distribution of gamma-rays (because the gamma-ray dis- 
tribution is additionally broadened due to deflections of 
electrons). The left panel of Fig. || shows SED of neu- 
trinos and antineutrinos received within different angles. 
The right panel of the figure presents the integral fluxes of 
neutrinos. The impact of the cutoff energy in the initial 
proton spectrum on the neutrino flux is demonstrated in 
Fig. g 

For comparison, the spectral energy distributions of 
protons, gamma rays and muon neutrinos are shown to- 
gether in Fig. |o| and Fig. |ll| for two distance to the 
source - 30 Mpc and 300 Mpc. 



IV. AN IMPULSIVE SOURCE: ARRIVAL TIME 
DISTRIBUTIONS 

Let assume that at the moment t = an impulsive 
spherically-symmetric source injects protons into the in- 
tergalactic medium. The multiple scattering of protons 
in the chaotic magnetic field results in the deviation of 
the motion of particles from the rectilinear propagation, 
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FIG. 6: Flux distributions of gamma rays observed within different angles in the direction of the source located at the distance 
r = 100 Mpc (left panel) and r = 300 Mpc (right panel). Calculations correspond to the IGMF B =1 nG and initial power- law 
distributions of protons with spectral index a = 2 and exponential cutoffs at Eo = 10 20 eV; 3 x 10 20 eV, and 10 21 eV. The total 
power of injection of protons into IGM is 10 44 erg/s. 



speed v p over the path S is 
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FIG. 7: Flux distributions of gamma rays observed within 
different angles in the direction of the cosmic ray source at 
the distance r — 100 Mpc. The calculations are performed for 
three different IGMF B =1 nG, 10 nG, 100 nG, assuming an 
initial power-law distribution of protons with spectral index 
a — 2 an exponential cutoff at Eq = 3 X 10 20 eV. The total 
power of injection of protons into IGM is 10 44 erg/s. The 
intergalactic absorption of gamma-rays due to interactions 
with EBL is not taken into account. 



therefore they arrive to the observer with significant time 
delays. The arrival time of the proton moving with a 



< = - = -+4.5xl0- 4 

v v c V 100 



10 1 



(57) 



For ultrarelativistic protons the second term is negligible, 
therefore in calculations we adopt v p = c. In this paper 
we will study the distribution of the arrival-time delays 
t = t — r/c ignoring the energy losses of particles. 

Let denote by P(t, r) d( dr the probability that the 
proton with arrival direction in the interval (£, £ + d() 
is detected at the distance r from the source in the time 
interval (r, r + dr). Here £ = 8 2 , where 9 is the angle 
between the proton direction at the point r and the vec- 
tor r. It is assumed that P s atisfi es to the condition 
of normalization given by Eq. (C27). The equation for 



the function P for a pulse of radiation in the small- angle 
approximation is obtained in Ref. [Q. In Appendix ^] 
we derive the exact relation between P and the standard 
distribution function /, and obtain P in a quite different 
(simpler) way than in Ref. j2J] . Namely, our treatment of 
the problem is based on the solution of equations written 
for the standard distribution function. 

Following to Ref. |24j we introduce the function G 
which is determined from the equation 



r 3 (6>2 



where the dimensionless parameters x and y are 
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FIG. 8: The fluxes of neutrinos observed within different angles in the direction of the cosmic ray proton source located at 
the distance r — 100 Mpc. The calculations are performed for the initial power-law distribution of protons with spectral index 
a = 2 and the exponential cutoff at Eo — 3 x 10 20 eV. The IGMF B — 1 nG, and the total power of injection of protons into 
IGM is 10 44 erg/s. Left panel is spectral energy distributions, right panel is integral fluxes. 




FIG. 9: The spectral energy distributions of muon neutrinos 
observed within different angles towards the source of cosmic 
ray protons located at r = 100 Mpc. The calculations are 
performed for initial power-law distribution of protons with 
a — 2 and three different values of the exponential cutoff: 
E = 10 20 eV, 3 x 10 20 eV, and 10 21 eV. 



Function G can be presented in the form of one- 
dimensional integral 



G(x 
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FIG. 10: The spectral energy distributions of gamma rays, 
muon neutrinos and protons observed within the polar angle 
3° from two identical source located at r = 30 Mpc (thick 
lines) and r = 300 Mpc (thin lines). The upper energy scale 
is for protons and neutrinos, the lower energy scale is for 
gamma-rays. The calculations are performed for the initial 
power-law distribution of protons with spectral index a = 2, 
the exponential cutoff Eo = 3 x 10 20 eV, and the total power 
of injection into IGM 10 44 erg/s. The IGMF is 1 nG. 
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FIG. 11: The spectral energy distribution of gamma rays, 
muon neutrinos and protons observed within polar angles 0.3° 
and 3° towards a source located at r = 30 Mpc. The param- 
eters for the IGMF and the proton spectrum are the same as 
in Fig. EH 



where z = \J s/(2i), j and j 1 are spherical Bessel func- 
tions: 

sinz . . sinz cosz ,„„. 
Jo («) = -— , 3l {z) = — 2 — . (62) 

The angular distribution of particles changes with 
time. It can be shown, by using Eqs. @ and @, that 

(9 2 )(r) =4cr/r, 

where 



(63) 



(6 2 )(r) = j6 2 G(x,y)dx/ jG(x,y)dx (64) 

o o 

is the mean square deflection angle at the moment r. 
Quite remarkably no model parameters enter in ( |63| ) in an 
explicit form. Thus, the measurements of 9 2 at different 
time periods allow an estimate of the distance to the 
source. This is a nice feature, because it could be the 
only channel of information about the distance to the 
source, if the latter is not active anymore. 

From Eqs. @ and @ follows that G(x,y) = at 
y < 0. We should note also the useful relation 



G(x,s) 



dy G(x, y) e 



(65) 



which allows us to obtain the moments of the function 
G: 



m n = I dyy n G(x,y) =i n —G(x,s) 



s=0 



(66) 



Let's write down the first three moments: 
m n = 3e~ 3x , 



771 , 



20 



(l + 2x) e 



-3x 



in 9 



2800 



(9 + 36a; + 28x 2 ) e 



-3.i' 



(67) 
(68) 

(69) 



Correspondingly the mean values for (y) and (y 2 ) are 



(y) = 



777,1 

m 



20 



(1 + 2x) 



2 7772 1 / 9 9 2 



(70) 
(71) 



For the dispersion of distribution A and the ratio -^p- we 
have 



&^(y 2 )-(y)' 



l 



1400 



(1 + 4x) , 



and 



A 

w 



2 1 + 4.t 2 
7 (1 + 2a;) 2 " 7 



(72) 



(73) 



This implies that we deal with a rather narrow distribu- 
tion. Rewriting Eq. (70) in the form 



c(r) = J_ 
20 



29 2 



(74) 



one can see that the measurement of (r) for the particles 
with different values of allows to estimate (0 2 ). 

Below we discuss two special cases of practical interest. 

A. Detection of protons with arbitrary arrival angles. 
This is the case discussed in Ref. p3]. In this case the 
distribution over t is described as 



oo 

/a- y^C,r)rfC = ^ ! yE(- 1 ) n " 1 ' 

n n—1 



with mean values for y: 



(y) = — , (v 2 ) = — , a = — 



(75) 
(76) 



B. Protons arriving along the radius- vector at the reg- 
istration point. For this case, substituting x = into 
Eq. (|60|), we obtain 



f B ^P(T,( = 0,r) = 



OO 9 



r 3 (6»2 



where < z\ < Z2 < ■ ■ ■ are the zeros of the function 
ji(z), located in the region z > 0. 

The functions f A and f R corresponding to Eqs. ( |75| ) 
and ( |77| ) are shown in Fig. |12]. 
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f A B , arb. units 




t , years 



FIG. 12: The arrival time distributions for the cases A (solid 
line) and B (dashed line) discussed in the text. The distance 
to the source is r = 10 Mpc, the energy of protons E = 
10 20 eV, and the strength of the magnetic field B = 1 nG. The 
curves are shown in arbitrary unites; for convenience they are 
normalized to 1 at the points of the maximum of distributions: 
max(/ AS ) = 1. 



As it follows from Eq. (pq), the arrival time r enters 
into P only in the form of combination of the variable y. 
Since 

lg y = lg r - 2 lg r - lg A - 2 lg B + 2 lg E + const , (78) 

the curves for other values of the relevant parameters, 
namely energy E, magnetic field B, correlation length A, 
and the distance to the source r, can be obtained by a 
simple shift along the r-axis. However, it should be noted 
that Eq. ([58]) is obtained in the approximation of ignor- 
ing the energy losses of protons. Therefore for the large 
distances, r > 100 Mpc, and especially for large energies, 
E > 10 20 eV, Eq. ( p8|) overestimates the arrival time, 
given that the energy of protons during their propaga- 
tion significantly exceeds the energy at the registration 
point (see Fig. 1). Therefore, for large distance Eq. ( pq ) 
should be treated as an upper limit for the time delay. 
On the other hand, since gamma-rays are produced at the 
very beginning of propagation of protons (within 10 Mpc 
or so), the curves calculated for a distance of order of 
10 Mpc, provide a quite accurate estimate for the arrival 
times of gamma-rays. 

V. SUMMARY 

In this paper the angular, spectral and time distri- 
butions of UHE protons and the associated secondary 
gamma-rays and neutrinos propagating through the in- 
tergalactic radiation and magnetic fields have been stud- 
ied based on the relevant solutions of the Boltzmann 
transport equation in the small-angle and continuous- 
energy-loss approximations. A general formalism for the 
treatment of the steady state distributions is provided 
in the form of relatively simple analytical presentations. 



The treatment of the secondary products, in particular 
the synchrotron gamma-radiation of electrons from pho- 
tomeson interactions is reduced to the consecutive appli- 
cation of the solutions which schematically can be pre- 
sented as 

Qp ^ fp ^ Qe ^ fe ^ Q~/ ^ f'y 

Here Qi denotes a source function and fi denotes a dis- 
tribution function. Q p is specified as spherically sym- 
metric source of protons. Q e is obtained from distri- 
bution function of protons as the final product of pho- 
tomeson interactions using the results Ref. [|l7j. Elec- 
trons generated in the pair production process of the 
first generation gamma-rays (from the decay of neutral 
7r- mesons) are also included in Q e . Finally, Q 7 corre- 
sponds to the synchrotron radiation of electrons with 
distribution function f e formed in the chaotic magnetic 
field. We consider the case of strong magnetic field, 
B > 10~ 9 G, when the electrons from photomeson in- 
teractions are cooled predominantly via synchrotron ra- 
diation. Such strong magnetic fields prevent the devel- 
opment of pair cascades at highest energies, and, at the 
same time, allow very effective conversion of the electro- 
magnetic energy released at photomeson interactions into 
synchrotron radiation. The latter peaks at GeV and TeV 
energies. The electromagnetic cascades are developed at 
lower energies at which the suppression of the Comp- 
ton cooling due to the Klein-Nishina effect is becoming 
more relaxed. These sub-cascades are initiated basically 
by the electrons-positron pairs produced at the inverse 
Bcthc-Hcitler process. However, because of deflections 
of low-energy electrons in chaotic IGMF, the gamma- 
rays produced during the cascade development lose the 
directionality. Moreover, if the initial energy distribu- 
tion of protons extends to 10 20 eV, the electromagnetic 
energy released in photomeson interactions greatly ex- 
ceeds the energy supply from the Bethe-Heitler process. 
On the other hand, the synchrotron radiation produced 
by highest energy secondary electrons not only provides 
an almost 100 % effective conversion into gamma-rays, 
but also preserves the initial direction of protons as long 
as the magnetic field does not exceed 10 -7 G. Remark- 
ably, while the main fraction of synchrotron gamma-rays 
and the highest energy neutrinos is produced in the prox- 
imity of the source, namely within the first ps 10 Mpc of 
the initial path of protons, the latter continue to suffer 
deflections with an enhanced rate (because of gradual de- 
crease of energy during the propagation through the 2.7 
CMBR), until they arrive to the observer. Therefore the 
gamma-ray and neutrino distributions appear to be more 
narrow than the angular distribution of protons. 

The distribution functions f p and f e are obtained by 
applying the Green function of transport equation to the 
source functions Q p and Q e , respectively. The angular 
part of f p and f e is a normal (Gaussian-like) distribution, 
the dispersion of which depends on the energy loss rate, 
the deflection angle per unit length and the distance to 
the source. / 7 is calculated by integration along optical 
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depth at different angles towards the source. 

For specific realizations of the scenario of small-angle 
deflection of charged particles, assuming that they move 
in a statistically isotropic and homogeneous turbulent 
magnetic field with Kolmogorov spectrum, we considered 
the IGMF in the interval from 10~ 9 to 10~ 7 Gauss and 
adopted 1 Mpc for the correlation length. The propa- 
gation of protons is considered, as long as it concerns 
the energy losses, as rectilinear with diffusion in angle. 
Transport of electrons is considered in the homogeneous 
magnetic field with random direction since their propa- 
gation length is of the order or less of 1 Mpc. 

Despite the small-angle scatterings, the related elon- 
gation of particle trajectories causes significant delays of 
their arrival time. The problem of propagation of parti- 
cles can be described by the steady state solutions if the 
lifetime of the source exceeds the delay times. Otherwise 
the problem should be treated as a time-dependent prop- 
agation of particles injected in the intergalactic medium 
by an "impulsive" source of extremely high energy pro- 
tons. This could be the case of solitary events like 
Gamma Ray Bursts or short periods (T < 10 5 year) 
of enhanced activity of active galactic nuclei. In this 
paper we discuss the case of an "impulsive" source, ig- 
noring the energy losses of protons. This approxima- 
tion limits the applicability of the derived time distribu- 
tion functions to the relatively nearby sources of protons 
located within 100 Mpc sphere of the nearby Universe. 
On the other hand, since the bulk of synchrotron radia- 
tion of secondary electrons is produced close the source, 
R < 10 Mpc, the time-dependent solutions derived for 
protons, can describe quite accurately the delayed ar- 
rival times of synchrotron photons from sources located 
at cosmological distances. 

The results presented in this paper for gamma-rays are 
valid for intergalactic magnetic fields in a specific (but 
perhaps the most realistic) range between B = 10 -9 - 
10~ 7 G. IMGF stronger than 10~ 7 G would lead to large 
deflections of charged particles, and thus violate the con- 
dition of small-angle approximation. On the other hand, 
IGMF weaker than 10 _10 G would reduce dramatically 
the efficiency of the synchrotron radiation since in this 
case the electrons are cooled predominantly via Compton 
scattering. The pair cascades initiated by these electrons 
also lead to GeV and TeV gamma-ray emission, how- 
ever these cascades form giant (hardly detectable) halos 
around the sources, unless the magnetic field is extremely 
weak, smaller than 10 _15 G. 

The realization of the scenario of synchrotron radia- 
tion of secondary electrons at the presence of a rela- 
tively modest magnetic field, B ~ 10 _9 G or larger, in 
the 10 Mpc proximity of the sources of highest energy 
cosmic rays, has higher chances to be detected, given 
the compact (almost point like) images at GeV and es- 
pecially TeV energies, and the very high (10 per cent 
or more) efficiency of conversion of the energy of pro- 
tons to high energy synchrotron gamma-rays. The fluxes 
of gamma-rays, protons and neutrinos shown in Figs. |7] 



- [n] are obtained assuming a power-law energy spec- 
trum of protons with a = 2 and total injection rate 
into IGM L p (> lGeV) = 10 44 erg/s. The expected 
gamma-ray fluxes are close to the sensitivities of Fermi 
LAT at GeV energies and the sensitivity of the Imag- 
ing Atmospheric Cherenkov telescope arrays at TeV en- 
ergies. While the total power of production of highest 
energy cosmic rays hardly can exceed, except for very 
powerful AGN, 10 44 erg/s, in the case of blazars with 
small beaming angles, the expected fluxes of gamma- 
rays could be significantly higher. Indeed, in the case of 
small deflections, the directions of injection of observed 
particles from the source are close to the observational 
line. Therefore the results for the spherically symmetric 
source remain valid also for the narrow jet with solid an- 
gle flbeam = , where D 2 > (6 2 ). Then the required 
power of the source (to be detected in gamma-rays) is 
reduced by a factor f2;, eam /47r which may be significantly 
small. 



Appendix A: The Green function for spherically 
symmetric source 

The Green function for spherically symmetric point 
source is obtained by integration of Eq. @ over all di- 
rections of the vector no. Let's rewrite Eq. (||) in the 
following form: 



n( j, j,, S(S(E,E )-r) ( A-Bn Q 

G(r,n,n ,E,E ) = _ T , m _, A — -exp' 



ce(E)ir 2 A 



(Al) 



Here 



A = 2{A 1 r 2 - 2A 2 r + A 3 ) - 2A 2 r{\ - n r n), 

B = 2{{A ir 2 - A 2 r)n r + (A 3 - A 2 r)n), (A2) 

where n r = r/r. Since the directions of n, no and n r 
are close to each other, 



Bn = \B\ cos e ~\B\ ( 1 - -j- 



(A3) 



Performing integration by saddle point method we ob- 
tain: 



exp 



A- Bn a \ 



2ttA 
TBI 



■ exp 



A-\B\ 



Taking into account 



n r n rjI , 

2 ' 



(A4) 
(A5) 



the expressions for A and \B\ can be written: 

A = 2Dr 2 + 9 2 A 2 r, 
\B\ = 2^{Dr 2 ) 2 - 9 2 (A ir 2 - A 2 r)(A 3 - A 2 r) (A6) 
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where Appendix B: Distribution function of electrons 

A 2 

D = A 1 -2-^ + -§. (A7) 

Expanding \B\ into series in terms of 9 to the second- 
order term in exponent and retaining the first term in 
denominator we find 

^<,M,*,-*§^£U(-£).<*, 

I 



After changing the order of integration in 

f e {r,n,E e ) = J Q e (f p (r ,n ,E p ))G(r - ro,n,no,E e ,E e o)dr Q no dE e0 (Bl) 

we arrive at the following integral over directions of the emission of electrons no a t the point ro and directions of r$: 

j f ( (n Q -n s ) 2 \ 
I= J°*»{ ~ J 

x exp(-^ l(r = r ° " |r " r °' no)2 " 2Mr r ° A |r " r °' no)(n " no) ± Mn no)2 ) dn n jn no , (B2) 

where 

r = rn r , r = r n s , r-r = r'n 07 (B3) 

A and D are defined in Eq. (0) and Eq. respectively. Taking into account that all directions are close, the 

integral can be presented in the following form: 

1 = j 'e-( A - B "-)dn„.dn„ , (B4) 

where 

A = X +X 1 (1- n r n ) + X 2 (1- n r n) +X 3 (1 - n n), 

B = Y\riQ + Y 2 n + Y 3 n r . (B5) 

Here we introduce the following notations: 



! + ^ ((r _ r ')2 +r 2 ); Xl = ^(A ir '-A 2 )r, X 2 = ^A 2 r, X 3 = |( 



Y 1 = ^ + ^(A 2 -A 1 r')r , Y 2 = -^A 2 r , Y 3 = ^A ir r Q . (B6) 
Since directions of n, no and n r are close, we can expand \B\ into series to the first order terms: 

\B\ ss v (1 — non) (1 — non r ) (1 — nn r ), (B7) 

V V V 

where 

I / = F 1 + y 2 + y 3 - (B8) 

The integration of / over dQ n3 by saddle point method gives 

fe-^-^dQ^ (B9) 
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To perform the integration over dQ no by the same method, we present the expression in the exponent in the following 
form: 

A - \B\ =Xo^+ (Xi + ™) (1 - n n r ) + (x 2 + (1 - nn r ) + (^X 3 + (1 - n n) 

Z 



Expanding \B\ into series 



where 



Z i Z2 Z$ 

Z Q + Z X + Z 3 + Z 2 (l - nn r ) - {Z x n r + Z 3 n) n Q = A - Bn . (BIO) 



\B\tt(i — (1-tmv), (Bll) 



H = Z x + Z 3 (B12) 



we find 



/« — e-C^-I^D. (B13) 

VfJ, 

Making replacements of all notations by their actual values and taking into account that 

r' « r - r , (B14) 

we finally obtain 

ir 2 AD ( flV \ 

1 = D e r 2 + D? o eXP {- D e r 2 + Drj ) ' (B15) 

where 

D e = Ax - 2-^ + ^| (B16) 

and is the angle between n r and n. The Integration over ro in the expression for f e (r,n,E e ) can be readily 
performed because of <5-function. 



Appendix C: Distribution of arrival times in the 
case of "impulsive" source 

In the case of spherical symmetry the proton distribu- 
tion function / = f(t, r, y) depends on time i, distance to 
the source r, and the variable fi — cos 9 = (nr)/r. Here 
n is a unit vector towards the direction of the proton 
speed. Let's normalize / using the condition 

oo 1 

Jdr J dnr 2 f(t,r,n) =1. (CI) 
o -1 

Then r 2 f(t,r, /i) dr dfj, is the probability that at the mo- 
ment t the proton is located in the layer (r, r + dr) and 
is moving in the direction within (/x, \x + dfi). Let as- 
sume that propagation of a single particle is fixed, i.e. 
the radius vector ro(t) and the direction no(t) are cer- 
tain functions of time. For this particle, the distributions 



over r and /j. are described by ^-functions: 

f (t, r,fi) = ^ S(r ~ r (t)) % - Mo (t)), (C2) 

where fx (t) = (n (t)r (t))/r Q (t). 

By averaging Eq. flC2| ) over the ensemble of particles 
gives the distribution function /: 

/(t,r^) = </o(t,r, M )). (C3) 

Let assume that for each particle ro(t) is a monoton- 
ically increasing function of time, i.e. there are no par- 
ticles in the ensemble with fi < 0. Then the equation 
r = ro(~t) has a unique solution with t — to(r), and thus 
one can write 

= ^ (< - <oW) - (C4) 
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Since in Eq. ( |C2| ) this expression is multiplied to 5(fi — 
fio(t)), in the denominator one can replace ^o(t) by /j, 
and take the factor 1/c/i out of the integral. This yields 



(<J(r-r (t))<J(M-W)(t))> 
= i<*(*-*o(r))*(/*-/io(r))), 

C/i 



(C5) 



where /2 (r) = Mo(*o( r ))- 

Function (5(t — to(r))5(n — Mo( r ))) has the meaning 
of the probability distribution for t and fi. Writing t = 
t + r/c, we obtain the probability distribution for r and 
fi at the point r: 

P(T,/*,r) = (*(T + r/c-t (r))*(M-/io(0))- (C6) 

From this equation follows that P satisfies the condition 
of normalization 



OO 1 

J dr j ' dfi P(t, r, fi) = 1 . 



(C7) 



Thus we arrive at the conclusion that the functions P 
and / are related as 

P{t, r ) = c ^r 2 f(r+r/c, r, n) = c^r 2 /'(r, r, fj,) . (C8) 

The distribution function satisfies the equation 

(C9) 



C OT 



where I is the collision integral. In the case of spherical 
symmetry 



(nV)/ = A* 



9/ , 1 - m 2 9/ 



9? 



(CIO) 



Replacing the variables in Eq. (C9) from (t, r) to (r = t — 
r/c, r) and presenting the collision integral in the Fokker- 
Planck approximation, we obtain 



l_A*5f + Of | 1 - /i 2 Of 
c dr dr r dfi 

= W±((i-f)°£\ 



(Gil) 



In the case of an impulsive source and no scattering 
(i.e. (0 2 s ) — 0) the distribution function normalized ac- 
cording to Eq. (pl|) is 



f(t, r, n) = 2ir S(r - ant). (C12) 
It is convenient to rewrite Eq. ( |C12| ) in the form 

f(t,r,n) = -^5(t-r/c)S(n-l). (C13) 

In order to demonstrate t hat the ge neral ized functions 
in the forms given by Eqs. ( |5l2|) and ( |cT^ ) are identical, 
one should multiply the right parts of these equations 



to an arbitrary function h(r,n) and integrate over the 
space coordinates and the direction of the vector n. This 
implies that at (9 2 ) =0 

/'( T ,r,/i) = -^*(r) S(jm-1). (C14) 

It is clear, from ge neral physical considerations, that in 
the limit r -> Eq. jSlj) is valid also at (6> 2 ) ^ 0. There- 
fore Eq. ( C14 ) can be treated as a boundary condition 

to Eq. (|C1H ) at the point r = 0. 

An analytical solution of Eq. ( pll ) is possible to de- 
rive in the small-angle approximation. In the case of 
multiple scattering, the average angle of deviation of at 
the distance r is of order of (r^ 2 )) 1 / 2 . Therefore for 
r -c l/(# 2 ) one can use the small angle approximation. 
Let fj, = 1 — C/2, and let us denote functio n /'(t , r, 1— £/2) 
by /'( T , r > 0- Assuming £ -C 1, from Eq. (Cll) we obtain 



To solve Eq. (G15) we apply the Fourier transformation: 



/W,C)= / /'(r,r,C)e 



dr. 



(C16) 



Function /' satisfies the equation 



df' iuC f/ 2(df> (f)2) dfdf\_ 



and the boundary condition given by Eq. ( |C14 ) becomes 



/'( w ,r,/i) = — 5(C), r- 
Let's search the solution in the form 



0. 



(C18) 



/' = e -<°M +6 < r > , (C19) 

where the fun ction s a(r), b (r) d o not depend on £. Sub- 
stituting Eq. ( p^| ) in Eq. ( p!7| ), we obtain the following 
ordinary differential equations: 



I ~0>- 



The solution to Eq. ( |C20[ ) is 

1 z i ( z ) 



a(r) 



r(0 2 ) 



(C20) 
(C21) 

(C22) 



where z = r us (8 2 ) / (2ic) . The arbitrary constant which 
appears in the solution is chosen requiring singularity at 
the point r = 0. At r -> the function a = 3/r(0 2 ). The 
solution to Eq. (C21) is 



6(r)=ln 3 



1 



r 3 j\(z) 



const . 



(C23) 
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therefore the function /' is defined as 



r 3 j^z) 



r(0 2 s ) 3i{z) 



(C24) 



For determination of the constant C on e should use the 
boundary condition given by Eq. ( |C18| ). In the limit of 
small r, and using the relation 

]i m (P e -W r )=6(Q, (C25) 



we find 



c r 3 ex Pl rm 



C^ S (Q. (C26) 



Comparing Eqs. ($32t ]and ( |C18|) , we obtain C = 2/c(6 2 s ) 
and then using Eq. (C8) we find P. In the small-angle 
approximation we can replace the factor /i in Eq. ( |C8|) by 
unity. In order to compare our results with the solution 
obtained in Ref. [pl| , we adopt C = l/c(9 2 ), which is 
equivalent to the change of the condition of normaliza- 
tion, namely instead of Eq. (C7) we use 



) oo 

dr [d(P(T,r,0 



(C27) 



where, because of rapid convergence, the upper limit of 
integration over d( is set infinity. In order to present the 
result in the form given by Eqs. (|8|) - (^), one should 
introduce, instead of the variable w, a new variable of 
integration s = ujr 2 (9 2 )/c. 



Appendix D: Emissivity function of synchrotron 
radiation in random magnetic fields 

For the case of chaotic magnetic fields one should av- 
erage out the standard formula for energy distribution of 
synchrotron radiation 



dN 1 _ V3 e 3 B (E^ 
dE^dt 2tt m e c 2 hE 1 \E C 



(Dl) 



where 



F(x) = x I K 5/3 (r)dr, E c = , (D2) 

' 2m„c 



over directions of magnetic field. After taking the per- 
pendicular to velocity component of magnetic field B± = 
Bsm9, where 9 is angle between B and v we come to 
the following double integral: 




FIG. 13: The emissivity functions for synchrotron radiation 
F(x) and G(x). The dashed line shows the ratio G(x) / F(x). 



After changing the order of the integration it can be writ- 
ten as a single integral 



G(x) = x I K 5/3 (t) 



(D4) 



that can be expressed in terms of modified Bessel func- 
tions: 

G{x) = ^r[(8 + 3x 2 )(n 1/3 ) 2 + xk 2/3 (2k 1/3 - 3xn 2/3 )}, 

(D5) 

where Ki/ 3 = K 1 / 3 (x/2), k 2 / 3 = K 2 / 3 (x/2). Note that 
while the function F(x) has a maximum at x = 0.2858 
(maxFfi) = 0.9180), the maximum of the function 
G(x) is shifted towards smaller values: x = 0.2292 
(maxG(x) = 0.7126). An alternative presentation for 
G(x) in terms of Whittaker's function has been derived 
in Ref. The functions F(x) and G{x), as well as the 
ratio G{x)/F{x) are shown in Fig. [i"3[ 

Although the function G(x) in Eq. (D5) is presented in 
a quite compact and elegant form, for practical purposes 
it is convenient to have approximation which does not 
contain special functions. Here we propose such approx- 
imations for F(x) and G(x) which provide an accuracy 
better than 0.2 % over the entire range of variable x: 



F(x) « 2.15x 1/3 (l + 3.06x) 1/6 
1 + 0.884a: 2 / 3 + 0.471a; 4 / 3 _ x 
1 + 1.64x 2 / 3 + 0.974a; 4 / 3 6 ' 



(D6) 



G(x) = [ sm0F( ^ = 1(f( -^-] sin 2 9d9. 
v ' J Vsin6>/ 4vr 2 J \sm9J 

o 

(D3) 



~. 1.808a; 1 / 3 1 + 2.21a; 2 / 3 + 0.347a; 4 / 3 

~ Vl + 3.4a; 2 /3 : 1 + 1.353a; 2 / 3 + 0.217a; 4 / 3 6 



(D7) 
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